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A large class of non-Markovian quantum processes in open systems can be formulated through 
time-local master equations which are not in Lindblad form. It is shown that such processes can 
be embedded in a Markovian dynamics which involves a time dependent Lindblad generator on an 
extended state space. If the state space of the open system is given by some Hilbert space H, the 
extended state space is the triple Hilbert space 7i ® which is obtained by combining the open 
system with a three state system. This embedding is used to derive an unraveling for non-Markovian 
time evolution by means of a stochastic process in the extended state space. The process is defined 
through a stochastic Schrodinger equation which generates genuine quantum trajectories for the 
state vector conditioned on a continuous monitoring of an environment. The construction leads to 
a continuous measurement interpretation for non-Markovian dynamics within the framework of the 
theory of quantum measurement. 
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I. INTRODUCTION 

An open quantum system is a certain distinguished 
quantum system which is coupled to another quantum 
system, its environment [ij. A particularly simple way 
of describing an open system is obtained in the Marko- 
vian approximation. In this approximation all memory 
effects due to system-environment correlations are ne- 
glected which usually leads to a Markovian master equa- 
tion, that is, to a linear first-order differential equation 
for the reduced density matrix p{t) of the open system 
with a time independent generator. Generally, one de- 
mands that the generator is in Lindblad form which 
follows from the requirements of the conservation of prob- 
ability and of the complete positivity of the dynamical 
map 0,a. 

A remarkable feature of Markovian master equations 
in Lindblad form is given by the fact that they allow 
a stochastic representation, also known as unraveHng, by 
means of a stochastic Schrodinger equation (SSE) for the 
state vector of the open system 0,IESH0|- A SSE gen- 
erates the time evolution of the state vector which results 
from a continuous monitoring of the environment of the 
system [mil. A specific realization {\ip{t)) ,t > 0} of 
the SSE is called a quantum trajectory: At each time 
t > the open system is known to be in a definite state 
\ipit)) under the condition that a specific readout of the 
monitoring of the system's environment is given. The 
reduced density matrix at time t is therefore obtained 
if one averages the quantity \'ip{t)) {^{t)\ over all possi- 
ble quantum trajectories. This means that the relation 
p{t) = E[\'ip{t)){il'it)\] holds, where the symbol E denotes 
the ensemble average or expectation value. 

In the Markovian case it is thus true that the envi- 
ronment acts as a quantum probe by which an indirect 
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continuous observation of the system is carried out. The 
description by means of a Markovian master equation 
in Lindblad form is, however, only an approximation 
which uses the assumption of short correlation times. 
For strong couplings and low temperature environments 
memory effects can lead to pronounced non-Markovian 
behavior. 

It is sometimes argued that the treatment of non- 
Markovian processes by means of master equations nec- 
essarily requires solving integro-differential equations for 
the reduced density matrix. Such equations arise, for ex- 
ample, in the application of the Nakajima-Zwanzig pro- 
jection operator technique 0, 0| which leads to dy- 
namic equations involving a retarded memory kernel and 
an integration over the past history of the system. 

However, the usage of another variant of the projection 
operator method allows in many cases the derivation of 
approximate or even exact non-Markovian master equa- 
tions for the reduced density matrix which are local in 
time. This method is known as the time-convolutionless 
(TCL) projection operator technique [il,|lE|lEli3- The 
non-Markovian character of the TCL master equation is 
refiected by the fact that its generator may depend ex- 
plicitly on time and may not be in Lindblad form. 

Time-local equations which are of the form of the TCL 
master equation have also been derived by other means, 
e. g. j by path integral and infiuence functional techniques 
]WL il^. A well known example is provided by the ex- 
act equation of motion for a damped harmonic oscillator 
coupled linearly to a bosonic reservoir [ioLElLE^. 

The fact that the TCL generator is generally not in 
Lindblad form leads to several important mathematical 
and physical consequences. In particular, a stochastic 
unraveling of the TCL master equation of the form indi- 
cated above does not exist: Any such process will auto- 
matically produce a master equation whose generator is 
in Lindblad form. The question is therefore as to whether 
one can develop a general method for the construction of 
stochastic Schrodinger equations for non-Markovian dy- 
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namics which do have a physical interpretation in terms 
of continuous measurements. It is the purpose of this 
paper to show that this is indeed possible. 

Our starting point is a time-local non-Markovian mas- 
ter equation for the density matrix p{t) on some Hilbert 
space Ti. with a time dependent and bounded generator. 
It will be demonstrated that the dynamics given by such 
a master equation can always be embedded in a Marko- 
vian dynamics on an appropriate extended state space. 
The non-Markovian dynamics thus appears as part of a 
Markovian evolution in a larger state space. 

If one chooses the extended state space as the Hilbert 
space of the total system, consisting of open system plus 
environment, this statement is of course trivial. How- 
ever, it turns out that the embedding can be realized in 
a fairly simple, much smaller state space, namely in the 
tensor product space TC C^. In physical terms this is 
the state space of a composite quantum system which 
results if one combines the original open system on the 
state space Tl with a further auxiliary three state sys- 
tem described by the state space C*^. The open system 
could be, for example, a damped quantum particle in- 
teracting with a dissipative environment. The auxiliary 
system can then be realized through an additional inter- 
nal degree of freedom of the particle which leads to a 
state space spanned by three basis states. 

It will be demonstrated that the dynamics in the ex- 
tended state space follows a Markovian master equa- 
tion with a time dependent generator in Lindblad form. 
The application of the standard unraveling of Marko- 
vian master equations to the dynamics in the extended 
state space therefore yields a stochastic unraveling for the 
non-Markovian dynamics. The resulting SSE generates 
genuine quantum trajectories which do admit a physical 
interpretation in terms of a continuous observation car- 
ried out on an environment. The construction thus gives 
rise to a consistent measurement interpretation for non- 
Markovian evolution in full agreement with the general 
setting of quantum measurement theory. 

The paper is structured as follows. Section ^ con- 
tains a brief review of the continuous measurement the- 
ory for Markovian dynamics. Time dependent genera- 
tors in Lindblad form are introduced in Sec. Ill AL and 
Sec. IIIBI treats the corresponding continuous measure- 
ment unraveling. The quantum measurement theory for 
non-Markovian evolution is developed in Sec. IIIII We 
introduce time-local non-Markovian master equations in 
Sec, nil A1 The embedding of these equations in a Marko- 
vian dynamics is constructed in Sec. IIII Bl whereas the 
derivation of the continuous measurement unraveling is 
given in Sec. IIII CI The construction of the SSE and 
its physical interpretation are illustrated by means of an 
example in Sec. IIVI 

A series of interesting stochastic unravelings of non- 
Markovian quantum dynamics is known in the literature. 
Section Fvl contains a discussion of our results and of the 
relations to alternative non-Markovian SSEs, as well as 
some conclusions. 



II. QUANTUM THEORY OF MARKOVIAN 
DYNAMICS 

A. Time dependent Lindblad generators 

We consider a density matrix W(t) on a state space H 
which obeys a master equation of the form: 



dt 



W{t) = C{t)W{t) 



^ -i[H(t),W{t)] + Y,Mt)w{t)jj{t) 

i 

-^E{4w^<w,w^w}- (1) 

i 

The commutator with the Hamiltonian H{t) represents 
the unitary part of the evolution and the Lindblad op- 
erators Ji{t) describe the various decay channels of the 
system. In analogy to the terminology used for classical 
master equations, the expressions JiWjj may be called 
gain terms, while the expressions {j/ji,T4^}, involving 
an anti-commutator, may be referred to as loss terms. 

Both the Hamiltonian H(t) and the operators Ji(t) are 
allowed to depend on time t. The generator C{t) of the 
master equation may thus be explicitly time dependent 
and does not necessarily lead to a semigroup. We ob- 
serve, however, that the superoperator Cito) is in Lind- 
blad form for each fixed to > 0. This means that £(to) 
is in the form of the generator of a quantum dynamical 
semigroup. The particular form of the generator derives 
from the requirements of complete positivity and of the 
conservation of the trace 0, Bl- 
under certain technical conditions which will be as- 
sumed to be satisfied here, one concludes that Eq. 
yields a 2-parameter family of completely positive and 
trace preserving maps V{t, s) (23,I3|- These maps can be 
defined with the help of the chronological time-ordering 
operator T as 



V{t,s) = Texp 



dTC{T) 



and satisfy 

V{t,s)V{s,t') = V(t,t'), 



t > s > 0, 



t> s>t' . 



(2) 



(3) 



In terms of these maps the solution of the master equa- 
tion at time t can be written as W{t) = Vit, s)W{s), 
where < > s > 0. Thus, V{t, s) propagates the density 
matrix at time s to the density matrix at time t. 

Each V(t, s) maps the space of density matrices into 
itself. This means that V(t, s) can be applied to any den- 
sity matrix W to yield another density matrix V{t, s)W . 
The range of definition of the maps V{t, s) is thus the 
space of all density matrices and is independent of time. 
Usually, one associates a Markovian master equation 
with a time independent generator. We slightly general- 
ize this notion and refer to Eq. Q as a Markovian master 
equation with a time dependent Lindblad generator C{t). 
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B. Stochastic unraveling and continuous 
measurement interpretation 

As mentioned already in the Introduction the master 
equation 1^ allows a stochastic unraveling through a ran- 
dom process in the state space 71. This means that one 
can construct a stochastic dynamics for the state vector 
in Ti which reproduces the density matrix W{t) 
with the help of the expectation value 



w{t)^E[\m){m\\- 



(4) 



To mathematically formulate this idea one writes a 
stochastic Schrodinger equation (SSE) for the state vec- 
tor . An appropriate SSE for which the expectation 
value Q leads to the master equation Q is given by 



d\^{t)) = -iG{\^{t)))dt 

jmm) 



E 



mt)\m)\\ 



dN,{t), (5) 



where we have introduced the nonlinear operator 



Gi\m)) 



\m)- (6) 



The term —iG{\^{t)))dt in Eq. ©, which is proportional 
to the time increment dt, expresses the drift of the pro- 
cess. This drift contribution obviously corresponds to the 
nonlinear Schrodinger-type equation 



dt 



(7) 



whose linear part involves the non-hermitian Hamilto- 
nian 



Hit) = Hit)-'-J2jjit)Mt) 



(8) 



The nonlinear part of the drift ensures that, although 
H{t) is non-hermitian, the norm is conserved under the 
deterministic time evolution given by Eq. (Q. 

The second term on the right-hand side of Eq. ^ rep- 
resents a jump process leading to discontinuous changes 
of the wave function, known as quantum jumps. These 
jumps are described here with the help of the Poisson 
increments dNi(t) which satisfy the relations: 



dNi{t)dNj{t) 
E[dN,it)] 



5,jdN,{t), 



'-dt. 



(9) 
(10) 



According to Eq. the increments dNi{t) are random 
numbers which take the possible values or 1. Moreover, 
if dNi{t) — for a particular z, we have dNjit) = for 
all j ^ i- The state vector then performs the jump 



m)) 



Mt)m)) 
mmm 



(11) 



Thus we see that Ni it) is an integer- valued process which 
counts the number of jumps of type i. 

We infer from Eq. ifTTHl that dNi{t) = 1 occurs with 
probability || J,(t)|$(t))|pd<. The jump described by 
ifTTll thus takes place at a rate of || Jj(t)|$(i))|p. The 
case dNi{t) = for all i is realized with probability 
1 — ||Ji(t)|$(i))|p(ii. In this case the state vector 
follows the deterministic drift described by Eq. 10. 

Summarizing, the dynamics described by Eq. Q yields 
a piecewise deterministic process, i. e., a random process 
whose reaHzations consist of deterministic evolution peri- 
ods interrupted by discontinuous quantum jumps. Since 
both the deterministic drift ((jjl as well as the jumps 
do not change the norm, the whole process preserves the 
norm of the state vector. 

The formulation of the dynamics by means of a SSE 
bears several numerical advantages over the integration 
of the corresponding density matrix equation Q. This 
fact was the original motivation for the development of 
stochastic wave function methods in atomic physics and 
quantum optics (for an example, see ^|). What is im- 
portant in our context is the fact that, additionally, the 
stochastic process given by the SSE allows a physical in- 
terpretation in terms of a continuous measurement which 
is carried out on an environment of the system. 

To explain this point we consider a microscopic model 
in which the open system is weakly coupled to a num- 
ber of independent reservoirs one reservoir for each 
value of the index i. Each reservoir Ri consists of bosonic 
modes hi\ which satisfy the commutation relations: 



t 1 



5 a 5, 



ijOXfj,- 



(12) 



The Lindblad operators Ji(t) appearing in the master 
equation couple linearly to the reservoir operators 



S.(i)=E 



■i{uJi - 



(13) 



where the uii are certain system frequencies, uJi\ is the 
frequency of the mode bi\ of reservoir Ri, and the gi\ are 
coupling constants. Thus, the Hamiltonian of our model 
is taken to be: 

Hj{t)=H{t) + ^J2 [Mt)Bl{t) + Jl{t)B,it)] . (14) 

We have included a factor I/VT, where F is a typical 
relaxation rate of the system which will be introduced 
below. The combination Ji/Vr is therefore dimension- 
less and the Bi have the dimension of an inverse time, 
choosing units such that h= 1. 

The time evolution operator over the time interval 
(t,t') will be denoted by U{t',t). The correlation func- 
tion of reservoir Ri can be expressed through the spec- 
tral density I{uj), which is assumed to be the same for all 
reservoirs: 



B,{t')Bl{t)\0) = / (ia;/(a;)e'('^' 



(15) 
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Here, |0) denotes the vacuum state defined by bi\\0) = 
for all i and all A. 

Suppose that the state of the combined system (open 
system plus reservoirs Ri) at some time t is given by 
\-qi{t)) = (g) |0). At time t' = t + t this state has 

evolved into the entangled state U {f , t)\'^ (t)) . We con- 
sider r to be a time increment which is small compared 
to the time scale of the systematic motion of the system, 
but large compared to the correlation time of the reser- 
voirs. Suppose further that at time t' a measurement of 
the quanta in the reservoir modes bix is carried out. Ac- 
cording to the standard theory of quantum measurement 
[2^ the detection of a quantum in mode bix projects the 
state vector onto the state 6-;^|0). The open system's 
state conditioned on this event thus becomes: 



1 



-.{0\b^xU{t',t)\^{t)), 



where 



(16) 



(17) 



is the corresponding probability. If, on the other hand, no 
quantum is detected, one has to project the state vector 
onto the vacuum state which yields the open system's 
conditioned state 



1 



where 



(o|t/(t',t)l*(t)), 



po^\mu{t',t)\^{tM^ 



(18) 



(19) 



is the probability that no quantum is detected. 

In the Born-Markov approximation the above expres- 
sion simplify considerably. We take a constant spectral 
density I{uj) = T/2tt, corresponding to the case of broad 
band reservoirs with arbitrarily small correlation times. 
We further neglect the Lamb shift contributions which 
lead to a renormalization of the system Hamiltonian. The 
expression ifTfiji then becomes (up to an irrelevant phase 
factor) : 



Mt)mt)) 
\Mt)\m)\\' 



(20) 



This expression is seen to be independent of A. Therefore, 
the total probability of observing a quantum in reservoir 
Ri is: 



Pi 



(21) 



which leads to the drift contribution —iG{\^{t)))dt of 
Eq. The probability for this event is found to be 



PQ 



= i-j2p.x = i-j2\\Mt)\m)\\'r. 



(23) 



Considering that the detected quanta are annihilated 
on measurement (quantum demolition measurement) we 
see that for both alternatives described above the con- 
ditional state vector of the combined system after time 
T is again a tensor product of an open system's state 
vector and the vacuum state of the environment. Thus, 
we may repeat the measurement process after each time 
increment r. In the limit of small r we then get a con- 
tinuous measurement of the environment and a resulting 
conditioned state vector of the open system that follows 
the SSE 

Summarizing, the SSE @ can be interpreted as result- 
ing from a continuous measurement of the quanta in the 
environment. This measurement is an indirect measure- 
ment in which the jump llll|l of the state vector |$(i)) 
describes the measurement back action on the open sys- 
tem's state conditioned on the detection of a quantum 
in reservoir Ri, while the nonlinear Schrodinger equa- 
tion yields the evolution of the state vector under 
the condition that no quantum is detected. The reaHza- 
tions of the process given by the SSE, i.e., the quantum 
trajectories thus allow a clear physical interpretation in 
accordance with the standard theory of quantum mea- 
surement. 



III. QUANTUM MEASUREMENT 
INTERPRETATION OF NON-MARKOVIAN 
DYNAMICS 

A. Time-local non-Markovian master equations 

We investigate master equations for the density matrix 
p{t) of an open system which are of the following general 
form: 



dt 



p{t) - mp{t) 



(24) 



+ [C^{t)p{t)Di{t) + D^{t)p{t)Cl{t)\ 

Oi 

-\ E {Di{t)Co.{t) + Cl{t)D^{t),p{t)] . 



The last two equation show that, conditioned on the de- 
tection of a quantum in reservoir Ri , the system state car- 
ries out the jump described in ljll|i . and that these jumps 
occur at a rate given by ||Ji(t)|$(i))|p [see Eq. ifTTHl ]. 

For the case that no quantum is detected expression 
III8II gives in the Born-Markov approximation: 



\m)-^Gmt)))r, 



(22) 



The Hamiltonian Hs{t), the Ca{t) and the Da{t) are 
given, possibly time dependent operators on the state 
space H. of the open system. The generator lC{t) may 
thus again depend explicitly on time. The master equa- 
tion II24II is however local in time since it does not con- 
tain a time integration over a memory kernel. The 
structure of JC{t) was taken to ensure that the hermitic- 
ity and the trace of p{t) are conserved. If we choose 
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Ca = Da = Ji/\/^ the generator /C(<) reduces to the 
form of a time dependent Lindblad generator C{t). The 
Markovian master equation is thus a special case of 
Eq. II24|I which will be referred to as non-Markovian mas- 
ter equation. 

With an appropriate choice for the Hamiltonian Hs [t) 
and the operators Ca{t) and Da{t), a large variety of 
physical phenomena can be described by master equa- 
tions of the form (I24II . For example, as mentioned in 
the Introduction a master equation of this form arises 
when applying the TCL projection operator technique 
[r3. IT^ to the dynamics of open system. The basic idea 
underlying this technique is to remove the memory kernel 
from the equations of motion by the introduction of the 
backward propagator. Under the condition of factorizing 
initial conditions one then finds a homogeneous master 
equation with a time-local generator JC(J,). The latter can 
be determined explicitly through a systematic perturba- 
tion expansion in terms of ordered cumulants 0, 0|. 
Specific examples are the TCL master equations describ- 
ing spin relaxation the s pin- boson model [2fll |. 
systems coupled to a spin bath [3nLl3ll |. charged particles 
interacting with the electromagnetic field ^21, and the 
atom laser [3^ . 

Moreover, several exact time-local master equations of 
the form II24II are known in the literature which have been 
derived by other means. Examples are provided by the 
master equations for non-Markovian quantum Brownian 
motion gy, [2l|, [22| and for the nonperturbative decay of 
atomic systems, which will be discussed in Sec. IIVI 

The existence of a homogeneous, time- local master 
equation requires, in general, that the initial state of the 
total system represents a tensor product state. For sim- 
plicity we restrict ourselves to this case since we intend 
to develop stochastic unravelings for pure states of the 
reduced open system. 

Due to its explicit time dependence the generator IC{t) 
of the master equation lf2ijl does of course not lead to a 
semigroup. But even for a fixed t > the superoperator 
IC(t) is, in general, not in Lindblad form, by contrast to 
the property of the generator £,{t) of the master equation 
QJ. To make this point more explicit we introduce op- 
erators Ea{t) — Ca(t) — Da{t) and rewrite Eq. II24|I as: 



d 
dt 



pit) - -i[Hsit),p{t)]+Viit)p{t)+V2it)p{t), (25) 



where we have defined the superoperators 



V,{t)p EE J2 



Cait)pciit)-^{ci{t)Ca{t),p} 



+Dait)pDi{t)-^{DUt)Dait),p} 
V,it)p EE -J2\Ec.it)pEiit)-l{Eiit)Eait\p} 



sign and violates therefore the complete positivity of the 
generator. 

Of course, we will assume in the following that the 
dynamics given by Eq. Ij24|l yields a dynamical map 
p(0) 1-^ p{t) for all times considered, that is, we sup- 
pose that Eq. (l2Hl describes the evolution of true density 
matrices at time t = into true density matrices at time 
t > 0. However, this assumption does not imply that 
the propagation of an arbitrary positive matrix at time 
t with the help of the master equation 1I2H1 necessarily 
leads to a positive matrix for future times. This is only 
guaranteed if we propagate a density matrix p{t) which 
results from the time evolution over the previous interval 
{0,t). 

A further important consequence of the form of the 
master equation Ij24ll is that it does not allow a stochas- 
tic unraveling of the type developed in Sec. rffRl Any 
unraveling of this kind would automatically produce a 
master equation with a time dependent Lindblad genera- 
tor. Trying to construct an unraveling which leads to the 
contribution 'D2{t)p{t) of the master equation, one would 
find a process with negative transitions rates, which is 
both unphysical and mathematically inconsistent. 



B. Markovian embedding of non-Markovian 
dynamics 

To overcome the difficulties in the development of a 
stochastic representation we are going to employ an in- 
teresting general feature of the non-Markovian master 
equation (I24II : Even if the generator K.{t) is not in Lind- 
blad form, it is always possible to construct an embedding 
of the non-Markovian dynamics in a Markovian evolution 
on a suitable extended state space. The precise formula- 
tion of this statement and its proof will be given in the 
following. 

The extended state space is obtained by combining the 
original open system on the state space Ti, with another 
auxiliary quantum system. The auxiliary system is a 
three state system whose state space C'^ is spanned by 
three basis states |2) and |3): 



span{|l),|2),|3)}. 



(26) 



The Hilbert space Ti, of the combined system then be- 
comes the triple Hilbert space 



(27) 



The extended state space is thus given by the tensor 
product of H and C'^, which in turn is isomorphic to the 
orthogonal sum of three copies Hi, TI2, Hs of Ti.. This 
means that states |$) in 7^ take the general form: 

1$) = 1^1) ® |l) + 1^2) ® |2) + I1/.3) ® |3), (28) 



One observes that I'i(t) is in Lindblad form, while (*) 
is not: The superoperator 'D2{t) carries an overall minus 



where IV'fe) € H for /c = 1, 2, 3. As a possible physical 
realization of the extended state space on may think of 
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Ti as the state space of a damped quantum particle with 
an additional internal degree of freedom which can be 
represented by a three level system. 

We now regard Eq. Q as an equation of motion on 
the triple Hilbert space, that is, W(t) is considered as a 
density matrix on Ti, governed by a master equation with 
time dependent Lindblad generator L(t). On the other 
hand, p{t) is a density matrix on TL satisfying the given 
non-Markovian master equation ll2Hl . We will assume 
in the following that the operators Ca{t) and Dait) are 
bounded. 



W{Q) = p(Q)(^\x){x\ 



dynamics in T-L 



W(t) 



dynamics in H. 



P(0) 



Pit) 



trWi2(t) 



Figure 1: Illustration of the embedding theorem: The initial 
density matrix p(0) evolves into p(t) according to the given 
non-Markovian master equation 12411 . This evolution can be 
embedded in a Markovian dynamics on the extended state 



space TC — Hi 



in which the density matrix 14^(0) evolves 



into W{t) following the master equation Q. 



where 



|X) = ^[|1) + |2) 



(32) 



is a state vector of the auxiliary system. We thus obtain 
W{0) by combining the open system in the state p(0) 
with the three state system in the pure state |x)(x|. It 
is obvious that W(0) given by Eq. Ij31|l is a true den- 
sity matrix on the triple Hilbert space H, i. e., we have 
W{0) > and trW{t) = 1. Moreover, Eq. ^ yields: 



Wi2{0) _ (1|T4^(0)|2) _ ip(0) 



trW^i2(0) tr(l|VK(0)|2) itrp(O) 



p(0), (33) 



which is Eq. (EHJ at time t — 0. 

To show that the relation II29II is valid for all times 
t>0 we have to demonstrate that the right-hand side of 
this relation satisfies the non-Markovian master equation 
II24|I . provided the Lindblad operators Ji{t) in Eq. (QJ 
are chosen appropriately. To simplify the presentation 
we first treat the case that the master equation 1I2H1 only 
involves a single operator C{t) and D{t). Thus we write, 
suppressing the time arguments. 



-p = -i[Hs,p] + CpD^ +DpC^ 



-Ud^c + c^d,p} 



(34) 



Our aim is to show that by an appropriate choice of the 
Hamiltonian H{t) and of the Lindblad operators Ji{t) in 
Eq. Q one can always achieve that the density matrix 
p{t) on Ti is connected to the density matrix W{t) on Ti. 
by means of the relation (tr denotes the trace): 



Pit) 



Wi2it) 

trWi2it)' 



Here, we have defined 



Wi2it) = {l\Wit)\2), 



(29) 



(30) 



which is an operator acting on H. We can regard this 
operator Wi2(i) as a matrix which is formed by the co- 
herences (off-diagonal elements) of Wit) between states 
from the subspace Hi and states from the subspace H2- 

This is the embedding theorem. It states that the non- 
Markovian dynamics of pit) can be expressed through the 
time evolution of a certain set of coherences W^it) of a 
density matrix Wit) on the extended space which follows 
a Markovian dynamics (see Fig. Q . 

To proof the embedding theorem we first demonstrate 
that the relation Ij29|l can be achieved to hold at time 
t = 0. If p(0) is any initial density matrix on H we define 
a corresponding density matrix on H by 



WiO) 



P(0)®^[|1>(1| 

p(o)® |x)(x|, 



|2)(2| + |1)(2| + |2)(1|] 




Figure 2: The extended state space U ® ^ TLi ® 7i2 ® H-i 
and the action of the operators Ji defined in Eqs. Ij35ll -ll38 | l. 

We define four time dependent Lindblad operators for 
the master equation (QJ: 



Ji(t) = Cit)®\l){l\+ Dit) ®\2){2\, 

J2it) = Dit)®\l){l\+Cit)®\2){2l 

hit) = r!(t)(^|3)(l|, 

Jiit) = r!(t)0 |3)(2|, 



(31) 



and the Hamiltonian: 

Hit) = Hsit)®[\\){\\ 
= Hsit)®h, 



|2>(2| + |3)(3|] 



(35) 
(36) 
(37) 
(38) 



(39) 
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where 1^ denotes the unit operator on the auxiliary space 
C'^. In Eqs. ll37ll and l|38|l we have introduced a time de- 
pendent operator Q.{t) on TL which will be defined below. 

According to the definitions Ii;-i5|l and the operators 
Ji and J2 leave invariant the subspace Tii © 7^2 which 
contains the states of the form iV'i) ® |1) + IV'2) ® |2). 
The operators J3 and J4 defined in Eqs. ll37jl and l(38jl 
induce transitions |1) |3) and |2) — > |3) between the 
states of the auxiliary three state system. The extended 
state space and the action of the Lindblad operators are 
illustrated in Fig. |2l 

We have from the above definitions (suppressing again 
the time arguments): 

= C^C®\l){l\ + D^D®\2){2\, (40) 

jIJ2 = D^D®\l){l\ + C^C®\2){2\, (41) 

Jg^Ja = |1)(1|, (42) 

4J4 = Vt^n®\2){2\, (43) 



and: 



4 

i=l 



4 J, = {n^n + C^C + D^D) ® (|1)(1| + |2)(2|) . 



(44) 

With the Lindblad operators defined in Eqs. ll35ll - l(38jl 
and with the Hamiltonian given by Eq. Ij39ll the master 
equation leads to: 



dt 



W,2 = -i{l\[H,W]\2) + Y,{l\J,Wjl\2) 



It- 

1=1 



l\{jjj,,W}\2). 



(45) 



On using Eq. Ij44|l we find 

4 

J2{M{JU^,w}\2) = {n^n + c^c + D^D,Wi2} 

i=l 

while Eqs. (1211)-® yield 

(l|JiW/|2) = CWi2D\ 
{1\J2W4\2) = DWi2C\ 



(IIJ3WII2) 



(l|J4Wi|2) =0, 



{l\[H,W]\2) - [Hs,Wi2] 



(46) 

(47) 
(48) 
(49) 
(50) 



Employing the relations lll6ll - lf5(11l in Eq. l(l5)l we arrive 



at: 



-i[Hs, W12] + CWuD^ + DWi2C^ 

-l{n^n + c^c + D^D,Wi2} ■ (51) 



Here we see the reason for our choice of the extended state 
space and of the Lindblad operators Ji. The operators 



Ji and J2 have been chosen in such a way that C acts 
from the left and from the right on the coherences 
W12, or vice versa [see Eqs. lUTjl and l(i8|l ]. On the other 
hand, since J3 and J4 induce transitions into the state 
|3) of the auxiliary system, the corresponding gain terms 
J3WJI and JiWj\ of the master equation (QJ do not 
contribute towards the equation of motion ll5T|l for the 
coherences W12 [see Eq. Ij49|l ]. The subspace 7^3 of the 
extended state space plays the role of a sink which will 
be used now to achieve that the loss terms of the master 
equation come out correctly. 

We observe that Eq. 115 l|l is already of a form which 
is similar to the desired master equation II34II . These 
equations differ, however, with respect to the structure 
of the loss terms which are given by the terms containing 
the anti-commutator. To get an equation of motion of 
the desired form we now choose SI to be a solution of the 
equation 

n'^n + C'^C + D'^D = aI + D'^C + C'^D, (52) 

which is equivalent to 

Vi^n = al - {C - D)\C - D). (53) 

Here, / denotes the unit operator on Ti. and a = a[t) is 
a time dependent non-negative number. Since Sl^fi is a 
positive operator a solution Vl of Eq. II53|I exists under 
the condition that the right-hand side of Eq. II53|I is also 
a positive operator. Thus, a must be chosen in such a 
way that 



a>\\{C-D)m? 



(54) 



for all normalized state vectors jV') in Ti. We note that 
it is at this point that the assumption of bounded oper- 
ators enters our construction. To make a definite choice 
we define a to be the largest eigenvalue of the positive 
operator (C - D)\C - D). This definition ensures that 
the inequality lf5ill is satisfied and that a solution Q of 
Eq. lf53jl exists. 

The solution of Eq. Ij53|l is, in general, not unique. 
If is a solution, then also Uft, where U is an arbi- 
trary unitary operator. Changing f2 into UV, does, how- 
ever, not infiuence the equation of motion ll5T|l since only 
the combination fi^Sl enters this equation. In the lan- 
guage of quantum measurement theory the transforma- 
tion p Slpf^'f is called a quantum operation j^^l- It 
describes the change of a density matrix p under a gener- 
alized measurement whose outcome occurs with probabil- 
ity tT{n^ilp). We can thus say that the unitary operator 
U, expressing our freedom in the choice of fi, affects the 
change of the system state p, but not the probability of 
its occurrence. 

Substituting Eq. lf52ll into Eq. ll5T)l we get: 

^Wi2 = -i[Hs,Wi2] + CWi2D^ +DWi2C^ 

-^{D^C + C^D,Wi2} -aWi2. (55) 
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We conclude from this equation that the trace of W12 
satisfies the equation 



— trWi2 — —atrWi2- 
at 



(56) 



Using this fact as well as Eq. one immediately 

demonstrates that the expression on the right-hand side 
of Eq. Ij29|l satisfies the desired master equation Ij34ll . 
which concludes the proof of the embedding theorem. 

The general case of an arbitrary number of operators 
Ca{t) and Da{t) in Eq. lf2ill can be treated in a similar 
way. To this end, one has to re-introduce the index a and 
to carry out a summation over a in the equations of mo- 
tion. Thus, for each value of a we have a corresponding 
ila{t) and an aa{t), as well as four Jia{t), i — 1,2,3,4. 

Finally we note that according to Eq. ll55ll the opera- 
tor W21 = {2\W\1) = satisfies the same differential 
equation as W12. Since also W2i(0) = Wi2(0), we con- 
clude that W2i{t) = Wi2{t) for all times. It follows that 
we can write for any operator A onH: 



tr{Ap{t)} = 



tr{{A^a^)W{t)} 
tT{{I^a,)W{t)} ■ 



where 



|1>(2| + |2>(1| 



(57) 



(58) 



is an operator on the auxiliary state space. This shows 
that the expectation value of all observables A in the 
state p(t) can be determined through measurements on 
the state W{t) of the extended system. 



C. Stochastic unraveling for non-Markovian 

processes 

The embedding of the previous section enables us to 
construct a stochastic unraveling for the non-Markovian 
dynamics given by the master equation lf2ijl . Since the 
master equation governing W{t) involves a time depen- 
dent Lindblad generator we can use the SSE developed 
in Sec. Ill Bl for this purpose. 

The SSE Q generates a stochastic process for the state 
vector !$(<)) in the triple Hilbert space. Employing the 
representation lf28jl we write 



\m) = 1^1 W> ® |1> + ® |2) + 1^3(0) ® |3). (59) 

As shown in Sec. IIIBI the density matrix W{t) on the 
extended state space is reproduced through the expecta- 
tion value W{t) = E[\<^{t)){^{t)\], and the norm of the 
state vector is exactly conserved during the stochastic 
evolution: 



mmt)) 

^ {Mt)\Mt)) 
= 1. 



In accordance with Eq. 113 l|l the initial state of the 
process is taken to be of the form: 

\m) = \^) ® ix) = ^ M ® ii> + 1^) ® |2)] , (61) 

where \ip) is a normalized random state vector in Ti, 
{tf\Lp) — 1, and Ix) is the fixed state vector of the aux- 
iliary three state system defined in l(32ll . We thus have 
iV'i(O)) = 1^2(0)) - -^W) and 1^3(0)) - 0. Hence, 
Eq. at time i = gives 

p{0)^EM{ip\]. (62) 

The embedding theorem now reveals that p{t) is ob- 
tained from the stochastic evolution with the help of the 
relation [see Eq. lf29jl ]: 

E{\Mt)){Mt)\) 



Pit) 



E{{Mt)\Mt))) 



(63) 



Thus we have constructed a stochastic unraveling for 
the non-Markovian dynamics. It is important to real- 
ize that our construction leads to an unraveling through 
a normalized stochastic state vector |$(i)) and that the 
process allows a definite physical interpretation in terms 
of a continuous measurement, as has been discussed in 

secnm 

We remark that for the case C = D the given mas- 
ter equation l(3ijl is already in time dependent Lindblad 
form. Our construction then yields a{t) = [see inequal- 
ity (jUl] and Ji = J2, as well as J3 — J4 = 0. This means 
that the jump operators Ji and J2 are identical and that 
the decay channels J3 and J4 are closed. It follows that 
l^sit)) = and that j-i/'i) evolves in exactly the same way 
as IV'2), that is, we have |V'i(0) = IV'2(0)- Equation l(63jl 
thus becomes p{t) = 2£'(|-0i(i)>(V'i(OI)- Note that the 
factor 2 is due to the normalization condition Il6()ll which 
yields (-01 j-i/'i) = (V'2|V'2) = 1/2- In the case C = D our 
construction therefore reduces automatically to the stan- 
dard unraveling of a master equation in Lindblad form. 



IV. EXAMPLE 

As an example we discuss a model for the non- 
Markovian decay of a two state system into the vacuum 
of a bosonic bath. The model serves to illustrate how 
to construct the embedding in a Markovian dynamics, 
and how to interpret physically the quantum trajectories 
generated by the resulting SSE. 



A. Construction of the process 

The interaction picture master equation of the model 
is given by 

'^-p{t) = -z[Hs{t),p{t)] (64) 



dt 



(60) 



+7(t) ( <J-p[t)a+ - i{a+a_,p(i)} 
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where 



(65) 



This is an exact master equation for the nonperturbative 
decay of a two state system with excited state |e) and 
ground state \g), which interacts with a bosonic bath 
cr+ ~ \e){g\ and ct_ — \g){e\ are the usual raising and 
lowering operators of the two state system. These oper- 
ators couple linearly to the bath through an interaction 
Hamiltonian of the form a-Q^{t) + a+Q{t), where Q{t) 
is a bath operator depending linearly on the annihilation 
operators of the bath modes. 

The real functions S{t) and j{t) are determined by the 
vacuum correlation function {0\Q{t)Q^ {ti)\0) of the bath. 
An example will be discussed below [see Eq. l|9Tll ]. The 
function S{t) describes a time dependent renormalization 
of the system Hamiltonian induced by the coupling to the 
bath (Lamb shift). Under the condition ^{t) > one can 
interpret the function j{t) as a time dependent decay rate 
of the excited state. But for certain spectral densities 7(t) 
may become negative in certain time intervals such that 
the master equation l(6ijl is not in Lindblad form and the 
generator is not completely positive. 

The master equation Ij64|l can however always be 
brought into the form of the non-Markovian master equa- 
tion l(3ijl by means of the definitions: 



C{t) 



m = 



To obtain the operator Q{t) introduced in Eq. 
first have to determine the quantity a = a{t) 
defined to be the largest eigenvalue of (C — 13)^ 
This operator is equal to for j{t) > 0, and 
2|7(i)|cr+CT_ for j{t) < 0. Thus we find 

{C-Dy{C-D) =aa+a- 

and 

a(i) = 17(01-7(0- 
Hence, Eq. Ij53|l takes the form 

= al - [C - D)\C - D) = a{I - 

which leads to an obvious solution: 

The Lindblad operators Ji defined in Eqs. 
therefore given explicitly by: 



Ji = [|l)(l|+sign7|2)(2|], 

J2 = [sign7|l)(l| + |2)(2|], 

.h = x/^(T+® |3)(1|, 

J4 = Vacr+ ® |3)(2|. 



(66) 

(67) 

lf53jl we 

which is 
{C-D). 
equal to 

(68) 



(69) 



(70) 

(71) 
l(38|l are 

(72) 

(73) 

(74) 
(75) 



B. Physical interpretation 

Considering times t for which ^{t) > 0, we have a = 
[see Eqs. l(69jl ] and, hence, J3 = J4 = 0. This means 
that for 7(t) > the decay channels described by J3 and 
J4 are closed: The process only involves the jumps of 
the state vector given by the operators Ji and J2. We 
infer from Eqs. (17211 and lf73jl that Ji and J2 induce down- 
ward transitions |e) — > \g) of the two state system (action 
of (T_). These transitions result from the projection of 
the system's state vector into the ground state \g) condi- 
tioned on the detection of a quantum in reservoir Ri or 

i?2. 

For 7(t) < the operators Ji and J2 again induce 
downward transitions of the two state system and, at the 
same time, introduce a relative phase factor of sign7 = 
— 1 between the states |1) and |2) of the auxiliary three 
state system. Moreover, in the case j{t) < the decay 
channels J3 and J4 are open: This enables the additional 
jumps of the state vector described by J3 and J4. Equa- 
tions iTtHi and ir75ll show that J3 and J4 lead to upward 
transitions \g) |e) of the two state system (action of 
tT+) with simultaneous transitions between the auxiliary 
states of the form |1) |3) or |2) |3). A re-population 
of the excited state |e) is thus possible through jumps into 
the auxiliary state |3), corresponding to the detection of 
a quantum in reservoir R3 or R4. 

A detailed analysis of the process can be given in terms 
of the statistics of the quantum jumps. To this end, we 
note that the waiting time distribution for the SSE © is 
given by 



Fih.to) - 1- 



Texp 



-I I dsH{s) 

'to 



, (76) 



where H is defined by Eq. 0. F(ti,io) is the probabil- 
ity that a jump takes place in the time interval {to,ti), 
given that the previous jump occurred at time to and 
yielded the state |$(io))- Note that this is a true cu- 
mulative probability distribution, i.e. F{ti,to) increases 
monotonically with ti and satisfies F{to, to) = 0. For our 
model we have [see Eqs. (|Hl, JSHl, (011 and lf52|l ]: 



H{t) = Hs{t)^l3 



(77) 



--(a(i)/ + 7(t)^+^-)®(|l)(l| + |2)(2|). 

Let us analyze the process starting from the initial 
state 



\m) 



® Ix), 



(78) 



and investigate the occupation probability Pg(t) of the 
ground state. From the master equation it is clear that 
this quantity is given by the simple expression: 



Pg{t) = 1 - CXP 



(79) 
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If 7(s) takes on positive and negative values, this is a non- 
monotonic function of time. Our aim is to illustrate how 
the stochastic dynamics reproduces this behavior and to 
explain the physical picture provided by the unraveling. 

In the stochastic representation we have the formula 
[see Eq. (l63|l ]: 



Pait) 



E[{g\Mt)){Mt)\9)] 
E[{Mt)\Mm ' 



(80) 



We denote the moment of the first jump by ti, the mo- 
ment of the second jump by t2. It follows from Eq. Ij76|l 
that the waiting time distribution for the first jump is 
given by 



F{ti,0) = 1 - exp 



ds\-f{s)\ 



(81) 



and that the waiting time distribution for the second 
jump, given that the first jump took place at time ti, 
becomes 



F{t2,h) = 1-exp 



dsa(s 



(82) 



We further denote the total number of jumps in the 
time interval {0,t) by N{t). Since the process starts 
from the state Ij78ll the first jump is given by the ap- 
plication of Ji or J2 which project the state vector into 
the ground state. Therefore, prior to the first jump we 
have 2{g\ipi{t)) {'ijj2{t)\g) — 0, immediately afterwards we 
get 2{g\'ipi{t)) {'ilj2{t)\g) — sign7(ii). During the time in- 
tervals in which 7(t) < a second jump described by 
J3 or J4 is possible by which the state vector leaves the 
manifold 1 7^2 and lands in 7^3 . Once the state vec- 
tor is in 7^3, no further jumps are possible. Of course, we 
get again 2(f;|^i(t)) (-02(0 Iff) = after the second jump. 
Summarizing we have three possible alternatives: 



N{t) = Q 2{g\Mt)){Mt)\9)^0, 

N{t) = l 2(.g|Vi(i)>(V'2(i)l5> =sign7(ti), 

N{t) = 2 ^ 2{g\Mt)){Mt)\9)^0. 

From these relations we find the expectation value 

E[2{g\Mt)){Mt)\9)] 



dtii^(ti,0)[l-F(t,ti)]sign7(<i) 



(83) 
(84) 
(85) 



(86) 



exp 



dtia{ti) 



exp 



dtih{ti)\ 



Here, the quantity dtiF{ti,0) is the probability that the 
first jump occurs in dti, while 1 — F{t,ti) is the proba- 
bility that no further jumps take place within {ti,t). 

A similar analysis can be performed to obtain the ex- 
pectation value of the quantity 2{'tp2{t)\tpi{t)) . One finds: 

N{t) = 2{MtMi(t)) = I, (87) 
Nit) - 1 ^ 2{Mt)\Mt)) - sign7(ti), (88) 
N{t)^2 2(V'2(i)IV^iW> =0. (89) 



Thus we get: 

E[2{Mt)\Mm 

= l-i^(t,0)+ / dtiF{ti,0)[l- Fit,ti)]sign-f{h) 
Jo 



exp 



dtia(ti) 



(90) 



The term 1 — F{t, 0) represents the no-jump probability, 
i. e., the contribution from the event N{t) = 0. The 
result l(9?Hl could have been obtained also directly from 
Eq. IHSl- Using finally Eqs. l|HS|l and ^ in Eq. JHOl we 
find, of course, the correct expression Ij79|l for the excited 
state probability. 




LU 
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40 



60 
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Figure 3: Simulation of the SSE © for the non-Markovian 
decay of a two-state system. Top: ground state probability Pg 
obtained from a sample of 10"" quantum trajectories (dots), 
and analytical solution (continuous line). Bottom: expecta- 
tion value l?[2('(/)2|-!/'i)]. Parameters: 70/A — 25, A/70 — 0.2. 



We illustrate the above analysis by means of an exam- 
ple. Figure 13 shows the results of a Monte Carlo simula- 
tion of the SSE @ for the case of a Lorentzian spectral 
density which is detuned from the transition frequency of 
the two state system by an amount A (damped Jaynes- 
Cummings model) . This leads to a bath correlation func- 
tion of the form 

{0\Q{t)QHh)\0) = l^e'Mt-t.)-^\t-t^\^ (gj) 

where is the Markovian relaxation time and is 
the correlation time of the bath. The simulation was car- 
ried out in a nonperturbative regime: While the Born- 
Markov approximation requires that 70/A ^ 1, the sim- 
ulation uses 7o/A = 25. In this regime j{t) becomes 
negative for certain time intervals. These intervals can 
be seen in the figure as those intervals over which Pg 
and i?[2(?/;2|V'i)] decrease monotonically with time. This 
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is a signature for the fact that transitions |1) |3) or 
1 2) — > 1 3) become possible through the channels J3 or 
J4. These channels are closed in the time intervals over 
which E[2{ip2\'4'i)] stays constant. 

The decrease of the ground state probability can be in- 
terpreted as due to virtual processes in which a quantum 
is emitted into the bath and re-absorbed at a later time. 
This is a clear non-Markovian feature of the dynamics. 
In the stochastic unraveling this decrease results from the 
contributions of those quantum trajectories which involve 
at least one jump and for which the first jump at time 
ti occurred during a phase in which j{ti) < 0. The first 
jump then yields a negative contribution to the expecta- 
tion values of 2{g\Tpi) {ip2\g) and of 2{^p2\ipi) as a result 
of the relative phase factor sign7(ii) — —1 between the 
states |1) and |2) introduced by Ji and J2 [see Eqs. ijMjl 
and Ij88ll ]. Moreover, in a possible second jump the state 
vector leaves the manifold Hi ® TI2 to end up in a state 
proportional to |e) (g) |3), which gives 2(^2 IV'i) = [see 
Eq. 

The decrease of Pg is therefore due to quantum tra- 
jectories for which a second jump is possible which leads 
to a re-excitation \g) — * |e). Thus we see that the vir- 
tual emission and re-absorption processes appear in the 
stochastic unraveling in the extended state space as cer- 
tain real processes, namely as jumps with Ji or J2 in- 
volving a negative phase factor (detection of a quantum 
in reservoir Ri or R2), and as jumps into the auxiliary 
state 1 3) (detection of a quantum in reservoir R3 or R4). 
This shows how the quantum memory effect of virtual 
emission and re-absorption processes is encoded and com- 
pletely stored in the continuous measurement record. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have developed a general method 
for the derivation of stochastic unravelings for non- 
Markovian quantum processes given by time-local master 
equations of the form 1I2H1 . The key point of the construc- 
tion is the fact that such master equations always allow 
a Markovian embedding in an extended state space with 
a rather simple structure, namely in the triple Hilbert 
space H'S)C^. Within this embedding the density matrix 
p{t) of the original open system is expressed through a 
certain set of coherences of the full density matrix W(t) 
on the extended state space. 

The transition to the extended state space can be 
viewed physically as the addition of a further degree of 
freedom which is realized by a three level system. This 
enables one to represent the given non-Markovian dy- 
namics by means of a suitable interaction with a Marko- 
vian environment consisting of the reservoirs Ri intro- 
duced in Sec. IIIBI Although the generator of the given 
non-Markovian master equation needs not be in Lindblad 
form, the corresponding dynamics in the extended state 
space is therefore governed by a time dependent Lindblad 
generator of the form of Eq. 1^ . The lifting to an appro- 



priate extended state space thus allows the derivation of 
stochastic Schrodinger equations for non-Markovian dy- 
namics through a consistent application of the standard 
theory of quantum measurement. The SSEs obtained in 
this way generate genuine quantum trajectories with the 
physical interpretation of continuous measurements. 

The construction of Sec. IIII B1 provides a fairly general 
method for the Markovian embedding of a given non- 
Markovian dynamics: Apart from the existence of the 
master equation and from the boundedness of the oper- 
ators Ca{t) and Da{t), no assumption was made regard- 
ing the interaction Hamiltonian, the spectral density, the 
reservoir state, its temperature, etc. The Markovian em- 
bedding could therefore be useful in itself since it enables 
one to employ well established and developed concepts 
from the theory of completely positive maps and Lind- 
blad generators in the study of non-Markovian master 
equations. 

Formulating a non-Markovian unraveling we made use 
of piecewise deterministic jump processes. In an electro- 
magnetic environment this corresponds, for example, to 
direct photodetection. It should be clear, however, that 
our derivation allows any unraveling in the extended state 
space. Alternatively one can use diffusion-type SSEs, 
which in a continuous measurement interpretation corre- 
spond to other detection schemes like homodyne or het- 
erodyne photodetection (icL ITll | . 

Various stochastic unravelings for non-Markovian dy- 
namics have been suggested in the literature, involving 
both jump processes [3^11^ as well as SSEs with col- 
ored noise ,37. .38. 39, 4Q|. The technique developed by 
Imamoglu is related to the method of pseudo modes 
[ill li^ l . It employs an approximate Markovian embed- 
ding of a given non-Markovian dynamics. This embed- 
ding is based on the assumption that the reservoir can 
be represented by means of an effective set of fictitious 
damped harmonic oscillator modes. The Markovian em- 
bedding of the present paper is realized in an entirely 
different way by the introduction of the triple Hilbert 
space, and avoids the expansion into pseudo modes. It 
should also be noted that in the present method the den- 
sity matrix p{t) is not given by the partial trace of the 
density matrix W{t) in the extended state space. 

A further interesting method has been formulated by 
Diosi, Gisin and Strunz (3^ l3^ . These authors em- 
ploy a nonlocal stochastic integro-differential equation 
for the state vector. As demonstrated by Gambetta and 
Wiseman it seems, however, that the nonlocal SSE 
does not admit a continuous measurement interpretation 
within the framework of standard quantum measurement 
theory (see also in this context). This means that 
measurements carried out at different times on the en- 
vironment will infiuence the dynamics in a way which is 
incompatible with the stochastic process. The SSE does 
therefore not generate genuine quantum trajectories in 
the sense it does for Markovian dynamics. 

A number of unravelings of non-Markovian dynamics 
has been proposed j4J,|4a| which are based on the idea of 
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propagating a pair \^pl{t)), \ip2(t)) of stochastic state vec- 
tors and of representing the reduced density matrix with 
the help of the expectation value p{t) = E[\il)i{t)) {ip2{t)\]. 
It is even possible to design an exact stochastic unravel- 
ing |3(l. l4q which neither requires the existence of a mas- 
ter equation nor a factorizing initial state. This method 
makes use of a pair of independently evolving product 
states in the state space of the total system. Related 
stochastic wave function methods have also been formu- 
lated for the description of bosonic and fermionic many- 
body systems 1471 Ml. and for the simulation of quan- 
tum gases (ifll l5nl | by use of the positive P-representation 
[51I l52 | . A measurement interpretation of these stochas- 
tic methods is however not available. 

A pair \^pi{t)), \^2{t)) of state vectors can be consid- 
ered as an element of the double Hilbert space 7i (8) = 
Ti-i © TI2 which is the tensor product of Ti. and the state 
space of a two state system. In Ref. 0| a stochastic un- 
raveling in the double Hilbert space has been constructed. 
Although this method has been demonstrated to provide 
a useful numerical tool, a continuous measurement inter- 
pretation seems again to be impossible. This is connected 
to the facts that not only the master equation in H but 
also the master equation in the double Hilbert space is 
generally not in Lindblad form and that the process does 
not preserve the norm of the state vector. 

We mention finally some restrictions of the present the- 



ory. Similar to the formulation of the Lindblad theorem, 
we made use of the assumption of boundedness of the 
operators Ca {t) and Da (t) in the non-Markovian master 
equation Ij24|l . This assumption excludes the immediate 
treatment of important cases, such as quantum Brown- 
ian motion which involves the unbounded operators for 
position and momentum of the particle. However, what 
is really needed in the proof is that the inequality Ij54|l 
is satisfied. Provided the dynamics of the state vector 
is confined to an effective subspace of H on which the 
right-hand side of inequality if^ill is bounded, we can still 
define a finite a(t) and construct the embedding. A fur- 
ther restriction of the theory is that for certain models 
time- local master equation of the form lf2ijl may not exist 
for very strong coupHngs. The latter can lead to singu- 
larities of the TCL generator and to a breakdown of the 
TCL expansion (an example is discussed in It is an 
important open problem whether a continuous measure- 
ment unraveling can be developed for such cases. 
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